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The history of star formation and mass assembly in early-type 
galaxies 
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ABSTRACT 

We define a volume limited sample of over 14,000 early-type galaxies (ETGs) selected from 
data release six of the Sloan Digital Sky Survey. The density of environment of each galaxy 
is robustly measured. By comparing narrow band spectral line indices with recent models of 
simple stellar populations (SSPs) we investigate trends in the star formation history as a func- 
tion of galaxy mass (velocity dispersion), density of environment and galactic radius. We find 
that age, metallicity and a-enhancement all increase with galaxy mass and that field ETGs 
are younger than their cluster counterparts by ~ 2 Gyr. We find negative radial metallicity 
gradients for all masses and environments, and positive radial age gradients for ETGs with 
velocity dispersion over 180kms _1 . Our results are qualitatively consistent with a relatively 
simple picture for ETG evolution in which the low-mass halos accreted by a proto-ETG con- 
tained not only gas but also a stellar population. This fossil population is preferentially found 
at large radii in massive ETGs because the stellar accretions were dissipationless. We estimate 
that the typical, massive ETG should have been assembled at z < 3.5. The process is similar 
in the cluster and the field but occurred earlier in dense environments. 

Key words: galaxies: elliptical and lenticular, cD, galaxies: formation, galaxies: abundances 



1 INTRODUCTION 

Observational determinations of the history of star formation in 
Early-Type Galaxies (hereafter ETGs) are of great importance be- 
cause hierarchical models of galaxy formation make firm predic- 
tions for the relation between age, metallicity and a-enhancement 
as a function of mass. These scaling relations, plotted against the 
observational proxy for mass, the velocity dispersion (hereafter a), 
have been the focus of many recent studies of ETGs (Annibali et 
al. 2007, Bernardi et al., 2006, de la Rosa et al. 2007, Gallazzi 
et al. 2006, Jimenez et al. 2007, Kuntschner et al. 2001/2, Lucey 
et al. 2007, Mateus et al. 2007, Nelan et al. 2005, Proctor et al. 
2004/8, Sanchez-Blazquez et al. 2006, Smith et al. 2007, Terlevich 
& Forbes 2002, Thomas et al. 2005). 

These studies clearly show that the simple picture of early for- 
mation of low mass galaxies, which then merge to form more mas- 
sive systems is incorrect. The oldest stellar populations are found 
in the most massive galaxies - one aspect of so-called "downsiz- 
ing". However, the observational constraint is sensitive to the epoch 
at which star formation ceased, not when it started, so that low 
mass ETGs can still be "old", as dynamically bound objects, but 
have a mean stellar age that is much younger. However, this is 
not sufficient to save the simple hierarchical picture because the 
a-enhancement is also seen to increase with mass - implying more 
rapid formation for massive objects. This is a prediction of mono- 



lithic collapse models. It is important to remember that both the 
star formation history and mass assembly history of ETGs deter- 
mine their evolution. 

The most common method of determining the age, metallic- 
ity and a-enhancement of ETGs is by comparison of the narrow 
band absorption line indices with simple stellar population (SSP) 
models. It is well-known that stellar population parameters based 
on these indices are also sensitive to minor episodes of recent star 
formation. Luminosity-weighted, SSP equivalent stellar population 
parameters, such as those discussed here, do not therefore, distin- 
guish between a genuinely young galaxy and an old galaxy that has 
experienced a "rejuvenation" event. 

Some recent studies, based on the colour-magnitude relation 
of ETGs using colours that are extremely sensitive to recent star 
formation, in fact paint a surprising picture of ETG evolution. 
Schawinski et al. (2007) have used GALEX ultraviolet imaging to 
show that 30% of massive ETGs show ongoing star formation and 
that this fraction is higher in low-density environments. A similar 
picture is given by mid-infrared Spitzer data. Both Clemens et al. 
(2008) and Bressan et al. (2006) find that ~ 30% of ETGs in the 
Coma and Virgo clusters have experienced some star formation in 
the recent past. 

A recent study by Rogers et al. (2007) has combined SDSS 
spectra and GALEX data to conclude that "weak episodes of re- 
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cent star formation" are a phenomenon more commonly associated 
with ETGs in the cluster environment, a result seemingly, but not 
necessarily, inconsistent with several studies that find older SSP- 
equivalent ages in denser environments. 

Here we repeat the analysis carried out in Clemens et al. 
(2006, hereafter Paper I), which used 3614 objects selected from 
data release 3 (DR3) of the SDSS. Applying the same selection 
criteria to DR6 we define a sample of 14353 ETGs, four times as 
many objects. 



2 SAMPLE SELECTION AND DATA ANALYSIS 

Sample selection is identical to that described in Paper I. The local 
environmental density is defined as the inverse of the distance to the 
fifth nearest neighbour, I/V5, corrected for the redshift dependent 
effect of survey boundaries. 

We measure the 21 line-strength indices of the original Lick- 
IDS system plus the additional indices HyF, H<5F, B4000 and HK. 
However, before measuring the narrow-band indices from each of 
the SDSS, spectra we smooth the spectra to the wavelength de- 
pendent resolution of the Lick-IDS spectra. This step is essential 
as the models that we use to derive the age, metallicity and a- 
enhancement are based on the Lick system. The index values are 
then corrected for the smoothing effects of the galaxy's velocity 
dispersion and aperture corrected to a standard normalized radius 
(a fraction of the half-light radius). See Paper I for a more detailed 
description. 

Besides being based on DR6, the present work differs from 
that of Paper I only in 2 respects. Firstly, in Paper I we chose to 
correct for the fixed angular diameter of 3" sampled by the SDSS 
fibre using the radial index measurements of 50 nearby E and SO 
galaxies (Rampazzo et al. 2005). Here, we make use of the r-band 
effective radii provided in the SDSS catalogue to derive the aperture 
correction directly. Statistically, a more reliable correction should 
be obtained in this way. However, this also means that we must 
choose a larger standard radius to which to correct (r e /10 was used 
in Paper I). Typically, the 3" diameter SDSS fibre samples ~ r e /2, 
with very few objects being so large that the fibre samples r e /10. 
To avoid extrapolating beyond measured radii we therefore aperture 
correct to a standard radius of either r e /2 or r e /4. 

A plot of index value versus the fraction of r e sampled by the 
fibre will show a gradient. However, the gradient is due not only to 
the radial gradients within each galaxy but also to the correlation 
between r e and a. Because index values are correlated with a this 
effect increases the magnitude of the gradient. To determine the 
aperture effect we therefore consider the variation of index value as 
a function of the fraction of r e sampled by the fibre in restricted bins 
of (T. In this way we minimize the effect of the index- o~ relations 
and determine the radial index gradients as a function of a. We plot 
index values as a function of r/r e in 5 bins, for 8 separate bins of 
C The bins are chosen to maintain a large number of objects in 
each bin and gradients are then derived by weighting each point by 
\j\fn where n is the number in the bin. The values of the radial 
index gradients and their variation with o~ are shown in Table Q] 

We use the radial index gradients to correct our index values 
for aperture effects. Firstly, we use the measured value of a for a 
given galaxy to determine the value of the radial index gradient Al : 

, dAl 

A/ = (O200 - 1) + A/ 200 (1) 

C/CT200 

where O200 ls me velocity dispersion in units of 200 kms~', and 



A/200 is the value of the index gradient for a = 200 kms . Values 
for dAI/d<J2oo and A/200 are given in TableQ] 

This value of the radial index gradient is then used with the 
measured value of r e of the galaxy to correct the index value to the 
equivalent radius, r = r e /2 or r e /4: 

/ e =/ + A/log(rr e /l"5) (2) 

where I c is the corrected index, / is the measured index value, r e 
is the effective radius of the galaxy in arcsec and r is the standard 
radius expressed as a fraction of r e (0.5 or 0.25 here). 

The use the the effective radii provided in SDSS (the Pet- 
rosian half-light radius) to effect the aperture correction has one 
caveat. That is, that these radii are not seeing corrected. Therefore 
for galaxies with a small angular diameter r e is over-estimated. Me- 
dian r-band seeing for SDSS imaging is l"4 (Adelman-McCarfhy 
et al. 2007, fig. 4) and the median effective diameter for our sam- 
ple is 5". We believe, however, that this is not a serious problem 
because the seeing also alters the light entering the spectroscopic 
fibre. The considerable agreement we find in values for the indices 
as a function of O with other studies (see below) reinforces this 
view. 

The second difference from Paper 1 is the removal of one 
index, G4300, from the fitting procedure used to derive the age, 
metallicity and a-enhancement as a function of o~ and environ- 
ment. This was done because the index seems less well modeled 
than previously thought. 



3 RESULTS 

Here we will use the index values to derive various evolutionary pa- 
rameters as a function of o~, environment and galaxy radius. Before 
that, however, we briefly evaluate various trends seen in the fully o~ 
and aperture corrected index values. 

3.1 Index values as a function of velocity dispersion 

In Table [2] we show the corrected index values as a function of 
<J. Because the values refer to indices aperture corrected to r e /2, 
they are not directly comparable to those of Paper I, where correc- 
tions were made to r e /W. Nonetheless, most indices show similar 
behaviour to those of Paper I. We briefly note here some of the 
larger differences seen in important indices (we refer to gradients 
expressed as dl /d\og(a) asa; and those as dlog//</log(o") asAi). 

C4668: The gradient of a\ = +4.4 much larger than that of 
Paper I (+1.8). This value is closer to the value of +5.2 found by 
Nelan et al. (2005, hereafter N05). Hj8 : The gradient of - 1 .0 (A t = 
—0.24) is shallower (—1.7 in Paper I). This value is in excellent 
agreement with Bernardi et al. (2003, hereafter B03) (Ai = —0.24) 
and N05 ( a x = -1.2). Fe5015: The gradient of a x = -1.5 in Pa- 
per I contrasts to the present value of + 1.6. This is much more con- 
sistent with N05 who find a\ = +1.0. Mg2: The value of a\ = 0.22 
is similar to that of Paper I, but is now more consistent with both 
B03 and Kuntschner et al. (2002). Mgb: The gradient of a x = 3.3 
( A; = 0.37) compares with 3.7 in Paper I. This is closer to that 
of B03 (A 1 = 0.32) and N05 (a x = 3.2). Fe5270: In Paper I a null 
gradient was found. The present value of a\ = +0.62 is consistent 
with N05 who find an identical value. 

We note, that globally, the new index gradients are much 
closer to those derived by N05 despite the fact that these authors 
aperture corrected index values to a fixed physical radius, rather 
than to a fixed fraction of r e as done here. 
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Index A/200 dM/daiaa S / N 



CNl[mag] 


-0.019 ±0.002 


-0.005 ± 0.008 


0.6 


CN2[mag] 


-0.018 ±0.002 


-0.010 ±0.009 


1.1 


Ca4227 


-0.113±0.014 


0.060 ±0.065 


0.9 


G4300 


-0.012 ±0.031 


0.379 ±0.149 


2.5 


Fe4383 


-0.423 ± 0.040 


0.600 ±0.188 


3.2 


Ca4455 


-0.121 ±0.017 


-0.008 ±0.081 


0.1 


Fe4531 


-0.188 ±0.030 


0.203 ±0.142 


1.4 


C4668 


-1.039 ±0.057 


-0.561 ±0.282 


2.0 




0.085 ±0.021 


-0.488 ±0.105 


4.7 


Fe5015 


-0.415 ±0.043 


-0.401 ±0.201 


2.0 


Mgl[mag] 


-0.024 ±0.001 


0.013 ±0.005 


2.5 


Mg2[mag] 


-0.028 ±0.001 


0.013 ±0.007 


2.0 


MgB 


-0.311 ±0.028 


0.286±0.128 


2.2 


Fe5270 


-0.173 ±0.023 


0.066 ±0.1 12 


0.6 


Fe5335 


-0.241 ±0.026 


-0.134±0.127 


1.1 


Fe5406 


-0.150±0.022 


0.010±0.108 


0.1 


Fe5709 


-0.053 ±0.015 


-0.052 ±0.074 


0.7 


Fe5782 


-0.103 ±0.015 


0.058±0.071 


0.8 


NaD 


-0.627 ± 0.037 


-0.546 ±0.173 


3.1 


Ti01[mag] 


-0.000 ±0.001 


0.004 ±0.002 


1.9 


Ti02[mag] 


-0.008 ±0.001 


0.002 ±0.002 


0.8 


B4000 


0.017 ±0.002 


-0.013 ±0.007 


1.9 


HK 


0.001 ±0.002 


-0.028 ±0.012 


2.4 


H<5F 


0.265 ±0.030 


-0.265 ±0.144 


1.8 


H7F 


0.290 ±0.031 


-0.558±0.154 


3.6 



Table 1. Radial index gradients as a function of velocity dispersion, a. The 
radial index gradients are expressed as, A/ = rf[og ^ r j . where / is the value 
of the index. The third column gives values for the dependence of the radial 
index gradients on o~ in units of 200 kms~', o"2oo- A/200 is the value of 
the index gradient at o~ = 200 kms~'. These parameters have been used 
in the aperture correction of all indices. Most indices show a well-defined 
radial gradient, but rather few show a convincing trend of this gradient on o~. 
The last column shows the ratio of the gradient in the third column and its 
eiTor and so is an estimate of the statistical significance of the gradient as a 
function of o". For those indices whose radial gradients show little trend with 
CT, the value of A/200 is a good measure of the radial gradient for galaxies 
of any a. 



3.2 Index values as a function of radius 

The spatial gradients we measure here describe the mean index 
value in apertures of varying radii, at fixed cr. This is in contrast 
with 'true' spatial index gradients, which are measured in increas- 
ing annuli projected on the galaxy. As a result, the values we mea- 
sure are smaller in magnitude than the true gradients. Our values, 
which we give in Table Q] are, however, directly applicable to aper- 
ture correction. 

All the narrow line indices, with the exception of HK, show 
a radial gradient. The gradients are negative except for the hydro- 
gen line indices and the 4000A break, B4000. Additionally, some 
indices show a well defined trend of index gradient with a (these 
can be quickly identified in the last column of Table[T|l. In all cases 
(except CN1 and CN2) the sense of this variation is that the in- 
dex gradient becomes less steep with increasing cr. In some cases, 
including Hj8, a significant index gradient at low values of a dis- 
appears completely for o~ > 250 kms -1 . We find no dependence on 
the radial index gradients with density of environment at fixed cr. 
The fact that some indices show gradients which decrease with in- 
creasing o~ does not necessarily imply that some process has acted 
to mix the stellar populations in more massive systems. We return 
to this below. 



Index 


a\ = dI/d\og(a) 


ao 


CN1 [mag] 


0.209 ±0.002 


-0.413 ±0.005 


CN2 [mag] 


0.224 ±0.003 


-0.416 ±0.006 


Ca4227 


0.400 ±0.020 


0.190 ±0.046 


G4300 


1.1 18 ±0.047 


2.758 ±0.106 


Fe4383 


1.759 ±0.061 


0.699 ±0.1 36 


Ca4455 


0.628 ± 0.024 


-0.172 ±0.053 


Fe4531 


1.1 84 ±0.044 


0.544 ±0.100 


C4668 


4.392 ±0.083 


-3.516±0.187 


H/3 


-1.033 ±0.032 


4.032 ±0.073 


Fe5015 


1.604 ±0.061 


1.242 ±0.1 37 


Mgl [mag] 


0.145 ±0.002 


-0.217 ±0.004 


Mg2 [mag] 


0.216±0.002 


-0.236 ±0.005 


Mgb 


3.305 ±0.057 


-3.325 ±0.128 


Fe5270 


0.617 ±0.038 


1.427 ±0.085 


rcjjjj 


n Q9? -i- n. rno 
u.yzz m u.uj7 


It U.Uo / 


Fe5406 


0.563 ±0.032 


0.384 ±0.072 


Fe5709 


-0.148±0.021 


1.206 ±0.047 


Fe5782 


0.348 ± 0.026 


-0.015 ±0.059 


NaD 


4.432 ±0.049 


-6.345 ±0.1 10 


TiOl [mag] 


0.029 ±0.001 


-0.032 ±0.002 


Ti02 [mag] 


0.045 ± 0.001 


-0.026 ±0.002 


B4000 


-0.128 ±0.003 


0.848 ± 0.006 


HK 


-0.064 ±0.003 


1.003 ±0.008 


H<5F 


-1.600 ±0.044 


4.032 ±0.099 


HyF 


-2.530 ±0.046 


4.333 ±0.104 



Table 2. Index values as a function of velocity dispersion, o~. The indices 
have been aperture corrected to a standard radius of r e j1. Straight line fit 
parameters are given as the gradient, a\ = d//rflog(o"), and intercept, ciq, 
of plots of index versus log o". 

3.3 Stellar population trends 

We now make use of the index values to derive the age, metallic- 
ity and a- enhancement of the galaxy population in our sample. 
We repeat the multiple linear regression procedure described in Pa- 
per I (to which the reader is referred for a detailed description). We 
briefly summarize the procedure here. 

Firstly, because our index values are not calibrated to the Lick 
system (due to the lack of Lick standard star spectra in SDSS) 
we consider index variations relative to the mean value at a cr of 
a = 200 kms~'. By working with these differential index values 
we avoid both the problem of absolute calibration to the Lick sys- 
tem and potential problems in the absolute calibration of the SSP 
models. We therefore perform a multiple linear regression accord- 
ing to equation 5 of Paper I. The linear regression is performed 
on the whole sample and on 2 subsets of environmental density, 
1/7-5 ^0.5 (typical of the 'field') and \/r$ ^ 1.5 (more typical of a 
cluster). We also perform the analysis on indices aperture corrected 
to 2 different radii, r e /2 and r e /4, to investigate radial trends within 
the individual galaxies. 

The results of this regression analysis showed that the car- 
bon abundance did not depend on a, having a constant offset as 
a function of environment, in contrast to Paper I. This difference 
is probably due to the better aperture correction used here and/or 
the exclusion of the G4300 index. We therefore remove the explicit 
carbon abundance from the regression analysis, allowing the carbon 
abundance to be included in the metallicity term. The simultaneous 
a posteriori fit of the model to three example indices is compared 
with the median of the data in the different bins of o~, in fig.Q] 

In the left panel of Figure [2] we show the results for the en- 
tire sample for two different radii. For r e /4 the trend of age with 
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a is very similar to that seen in Paper I with a steady rise in age 
from the lowest mass systems and an approximately constant age 
for galaxies with o~ > 230 kms~'. For the larger aperture, how- 
ever, the trend is slightly different, with a less pronounced flatten- 
ing towards high values of a. For o~ > 300 kms -1 the mean age 
is ~ 0.05 dex older for r e /2 compared to r e /4. This corresponds 
to an age difference of ~ 1 Gyr for a galaxy of age 10 Gyr. The 
cross-over point of the 2 lines in the top panel of fig. [2] shows that 
galaxies with o~ > 180 kms _I have positive radial age gradients. 
Neither Sanchez-Blazquez et al. (2007) nor Mehlert et al. (2003) 
find evidence of radial age gradients in ETGs. We derive an age 
gradient in massive galaxies despite the absence of a gradient of 
the Hp index. 

The increase of metallicity with o~ is less strong than that 
found in Paper I. The gradient for the indices, corrected to r e /4, for 
o~ > 160 kms -1 is d log(Z)/d log(c>) ~ 0.42. This value is simi- 
lar to that found by Nelan et al. (2005), Thomas et al. (2005) and 
Smith et al. (2007) but smaller than that of Kuntschner et al. (2001) 
and Graves et al. (2007). For o~ < 160 kms _I , however, there is 
no significant trend of index value with a. For the r e /2 aperture, 
the metallicity is ~ 0.05 dex lower so that early-type galaxies are 
less metal rich at larger radii. Negative metallicity gradients have 
also been reported by Proctor et al. (2008), Annibali et al. (2007), 
Sanchez-Blazquez et al. (2007) and Harris & Harris (2002). There 
is also evidence that the metallicity gradients are steeper for more 
massive galaxies as found by Forbes et al. (2005). The metallic- 
ity gradient likely compensates the age gradient to remove radial 
gradients in indices like H/3. 

The trend of a-enhancement with o~ is also slightly less steep 
than found in Paper I, with d [a/Fe]/d log(o") ~ 0.55, similar to 
that found by Annibali et al. (2007) but steeper than the ~ 0.3 found 
by several authors ( Thomas et al., 2005, Kuntschner et al., 2001, 
N05, Bernardi et al., 2006, Smith et al., 2007). Although most sim- 
ilar studies refer to a smaller radius, we see that the larger aperture, 
r e /2, has a marginally shallower gradient. The a-enhancement 
within this larger aperture is slightly lower, with the difference 
being largest (0.03 dex) for the most massive galaxies. This does 
not support the 'outside-in' ETG formation scenario (Pipino, Mat- 
teucci, & Chiappini, 2006). A negative a-enhancement gradient is 
also seen in the halo stars of the Galaxy (Fulbright, 2000). 

In the right panel of Figure [2] we show the variation of evo- 
lutionary parameters as a function of the density of environment 
for index values aperture corrected to r e /4. In both environments 
the increase in age with a is similar, with a flattening above 
o" ~ 200 kms~' . objects in dense environments (l/r$ ^ 1.5), how- 
ever, are 0.087 dex older than those in less dense environments 
(I/V5 ^ 0.5). This is a difference of 2 Gyr if the ages are close 
to 10 Gyr and is consistent with several earlier studies (Terlevich 
& Forbes 2002, Kuntschner et al. 2002, de la Rosa et al. 2007, 
Sanchez-Blazquez et al. 2006). The flattening of the age-o" relation 
is also slightly more pronounced for the lower density environment. 
Similar age trends were seen in Paper I. 

There is marginal evidence that the metallicity is lower in high 
density environments. Formally the difference is 0.020±0.013 dex. 
Thomas et al. (2005) also found a small environmental dependence 
on the metallicity in the same sense and de la Rosa (2007) finds 
a difference of 0.11 dex between Hickson Compact Groups and 
the field. Other authors have found both larger differences in the 
same sense (Proctor et al. 2004, Kuntschner et al. 2002), no effect 
(Bernardi et al., 2006, Annibali et al. 2007) and the opposite effect 
(Gallazzi et al. 2006, Mateus et al. 2007). In Paper I no difference 
was found. 
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Figure 1. A posteriori comparison of the models with selected indices (only 
three are shown but all indices are similarly reproduced). The thick red lines 
are the simultaneous solution to the whole set of indices (not fits for each 
index), while the blue lines trace the median data value. 
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Figure 2. Age, metallicity and a-enhancement variations as a function of 
O". A linear regression analysis has been carried out simultaneously on the 
H/3, H<5, Mgl, Mg2, MgB, Fe4383, Fe4531, Fe5270, Fe5335 and C4668 
indices. Left: Radial variation. The two lines in each plot refer to index 
values corrected to r e /4 (solid line) and r e /2 (dashed line). Right: Effect 
of environment. The solid line represents the entire sample, diamonds only 
those objects in low density environments (l/r^ < 0.5) and triangles only 
those in high density environments ( 1 / rs > 1 .5). Values are differences with 
respect to those of the entire sample at a = 200 kms~' . The central o" bins 
typically contain > 10 3 galaxies, whereas the highest bin contains only 16. 



No environmental effect is found for the a-enhancement, in 
agreement with Kuntschner et al. (2002), Thomas et al. (2005), An- 
nibali et al. (2007) and Gallazzi et al. (2006). However, Proctor et 
al. (2004), Bernardi et al. (2006) and Lucey et al. (2007) all find 
increased a-enhancement in denser environments. 



4 DISCUSSION AND CONCLUSIONS 

We find positive correlations between age, metallicity and a- 
enhancement and the velocity dispersion, o", in ETGs. Galaxies in 
dense environments are ~ 20% older than those in low density en- 
vironments for all o~ (~ 2 Gyr for an age of 10 Gyr). The trend 
with age flattens above o~ ~ 200 kms~' , especially for galaxies in 
low density environments. We find a marginally significant trend 
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towards higher metallicities in low density environments but the 
environment has no effect on the a-enhancement. 

Apart from the marginal metallicity difference between field 
and cluster the results are very similar to those of Paper I. There 
we concluded that an anti-hierarchical scenario, in which star for- 
mation lasts longer but with lower efficiency in lower mass objects 
(see Granato et al. 2004) was consistent with the data. Here the ad- 
ditional determination of SSP parameters as a function of galactic 
radius places additional constraints on the evolutionary scenario. 

Massive ETGs (a > 180 kms -1 ) have positive radial age gra- 
dients, negative metallicity gradients and marginally significant 
negative a-enhancement gradients. When a massive halo becomes 
non-linear it accretes smaller halos which started to collapse at ear- 
lier times. The radial trends suggest that these halos do not contain 
only gas, but also pristine stars. The gaseous component falls dis- 
sipatively into the potential well of the massive (proto-)spheroid, 
fueling rapid star formation. The increase in mass increases the rate 
and efficiency of star formation, driving the main correlations with 
galaxy mass. The pristine stellar component of each sub-halo, how- 
ever, being dissipationless, is deposited at a radius consistent with 
the angular momentum of the encounter. These stars, which are 
slightly older, more metal poor and have moderate a-enhancement 
will therefore be spread over larger radii. At early times, rapid gas- 
rich mergers lead to an almost monolithic formation, at later times 
mergers become increasingly "dry". Very similar scenarios have 
been proposed to explain both the bi- modal metallicity distribution 
of globular clusters and the greater radial scale length of metal- 
poor relative to metal-rich globular clusters in elliptical galaxies 
(Cote, Marzke & West, 1998, Bekki et al. 2008). Our results imply 
that the metal-poor globular cluster population in ETGs should be 
older, less metal-rich and slightly less a-enhanced than the metal- 
rich clusters. 

Because the age difference at larger radii is actually the lumi- 
nosity weighted SSP equivalent age in a larger aperture (not an an- 
nulus), the value of 0.05 dex, ~ 1 Gyr, is a lower limit to the real age 
difference at larger radii. This time difference limits the assembly 
redshift of massive ETGs simply due to the lack of time to accom- 
modate the formation of stars in the lower mass halos. Our limit 
translates into an upper limit to the assembly redshift of massive 
ETGs of z < 3.5; in which case the stars in low mass halos formed 
at z ~ 10, for a standard cosmology (Hq = 70, Q. m = 0.3, £2a = 0.7). 
This also provides an estimate of the star formation rate in the as- 
sembled spheroid. If the final stellar mass were 10 12 M then stars 
must have formed at a rate ~ 10 12 M / 1 Gyr ~ 10 3 M© yr~ 1 . 

We conclude by stressing the statistical nature of our results. 
Because a galaxy's velocity dispersion is a function of both the halo 
mass and virialization redshift, variations in these parameters may 
render small samples insensitive to the trends we find. 

The catalogue on which this article is based can be found at, 
www .mrao . cam. ac . uk/~bn204/galevol/ clemensetal08 .html. 
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